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Correlations and the formation of bound states (nuclei) are essential for the properties of nuclear 
matter in equilibrium as well as in nonequilibrium. In a quantum statistical approach, quasipar- 
ticle energies are obtained for the light elements that reflect the influence of the medium. We 
\Q ■ present analytical fits for the quasiparticle energy shifts of light nuclei that can be used in various 

Cn| \ applications. This is a prerequisite for the investigation of warm and dense matter that repro- 

duces the nuclear statistical equilibrium and virial expansions in the low-density limit as well as 
relativistic mean field and Brueckner Hartree-Fock approaches near saturation density. 
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1. Introduction 

(N 

We consider nuclear matter at moderate temperatures, T < 20 MeV, and subsaturation den- 
sities, i. e. baryon densities n < n s sa 0.16 fm . The asymmetry given by the proton fraction 
\ Y p = n p /n is arbitrary, where n p denotes the total density of protons, and n n -n-n p v& the total 

t^J" ' density of neutrons. In this regime, correlations between nucleons are important, in particular 

the formation of bound states, i.e. nuclei. 

A general treatment of correlations and cluster formation of the interacting many-nucleon 
system can be given within a quantum statistical framework [ 1]. Equilibrium properties such as 
the equation of state (EOS) as well as transport properties are fundamentally influenced by the 
composition of nuclear matter. The investigation of the properties of nuclear matter in this region 
of parameter values for T, n, and Y p (warm nuclear matter) is of relevance for astrophysical 
applications such as supernova explosions, but also for the description of heavy ion collisions 



S ' ( fflC )- 



Near saturation density, a quasiparticle approach can be used to describe nuclear matter. 
The medium effects result in a self-energy shift of the nucleons that can be incorporated in 
the chemical potentials fi„,fi p of neutrons and protons, respectively, and in the effective masses 
m* n ,m p . Different semi-empirical approaches such as the Skyrme parametrization ]2[] or more 
sophisticated relativistic mean-field (RMF) approaches J^t] are used to describe the quasiparti- 
cle energy shift of nucleons in dense matter, in particular in astrophysical applications. A more 
fundamental treatment is possible in quantum statistical approaches such as the Dirac Brueckner 
Hartree-Fock (DBHF) approximation, see [4]. Medium effects become relevant at nucleon den- 
sity above ~ 10~ 4 frrT 3 . However, such single-nucleon quasiparticle approaches fail to describe 
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warm nuclear matter at low densities because correlations, in particular bound state formation, 
are neglected in mean-field theories. 

In the low-density limit, warm nuclear matter can be described as a mixture of different com- 
ponents, the clusters. In equilibrium, a mass action law can be derived, and the composition is 
determined by nuclear statistical equilibrium (NSE), see [5]. The full account of two-particle cor- 
relations (A = 2) in the equation of state results in the Beth-Uhlenbeck expression for the second 
virial coefficient [6]. However, NSE in which the interaction between the different components 
is neglected, is possible only below nucleon densities of ~ 10~ 4 fm~ 3 . To include interactions 
on the mean-field level, a generalized Beth-Uhlenbeck formula has been derived in a quantum 
statistical approach flO] that reproduces the exact low-density limit as well as the quasiparticle 
approach at saturation density. In addition to the single-nucleon self-energy shifts, also Pauli 
blocking has to be considered to describe the mean-field effects for bound states consistently. 

The formation of higher clusters (A > 2) can also be included in the quantum statistical 
approach [8i]. A cluster-virial expansion H leads to the generalization of the Beth-Uhlenbeck 
formula for interacting clusters that are considered as new components (chemical picture). In 
particular, the formation of a particles is of importance. This has been accounted for in astro- 
physical equations of state |2l[3|]. The contribution of clusters near saturation density, however, 
has to be suppressed. In a simple phenomenological approach, this suppression is modeled in 
an ad hoc way by an excluded volume picture. A more systematic treatment can be given in a 
quantum statistical approach where single-nucleon self-energy shifts as well as Pauli blocking 
are included . 

Recently, interest arose in a systematic treatment of the nuclear matter EOS [ 10]. A combi- 
nation of NSE with the RMF approach was given in Ref. Ill ill treating the medium effects for 
clusters via the excluded volume picture. Based on a quantum statistical approach, the nuclear 
matter EOS including light cluster (A < 4) formation was given in Ref. 01211 . The mass fraction 
of different clusters has been calculated, and it was shown how clusters are suppressed with in- 
creasing density due to Pauli blocking. This quantum statistical approach that reproduces NSE 
and virial expansions in the low-density limit has been applied to supernova explosions lfl3ll as 
well as to HIC to determine the symmetry energy in the low-density region II 1411 . 

In the quantum statistical approach, the interaction between the different clusters is treated 
by the method of thermodynamic Green functions HI. Not only the free nucleons are treated 
as quasiparticles with momentum dependent energy shifts, but also the clusters are described as 
quasiparticles with energy shifts depending on the center of mass momentum P of the cluster. 
The composition of warm nuclear matter as well as the macroscopic properties are determined 
by the cluster quasiparticle energy shifts. The aim of the present work is to provide the medium 
modification of the energies of light elements (A < 4) in dense nuclear matter in a compact form, 
that can be easily used in further calculations. We generalize and improve previous parametriza- 
tions of the quasiparticle energies of light elements B15I1 as a function of P, T, n, and Y p . Our goal 
is to achieve an accuracy of better than 1% in a wide region of parameter values. 

The present work gives the quasiparticle energy shifts only for the light elements. Heavier 
elements need further investigation and are subject to future work. They are of importance at low 
temperatures. Regions of the phase diagram that are dominated by light elements are shown, e.g., 
in ll ill . We also neglect weak interaction effects that lead to j3 equilibrium, as well as Coulomb 
interaction. A homogeneous background of electrons is assumed to compensate the charge of 
protons. The homogeneous nuclear matter considered here may become thermodynamically 
unstable against phase separation. Then, large droplets and more complex structures can be 
formed, and Coulomb interaction has to be taken into account (see M 1 611 and references therein) 
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which, however, is not included in our treatment. Also the formation of quantum condensates 
(pairing, quartetting) will not be discussed here. 



2. Quantum statistical approach to warm nuclear matter 

We shortly sketch the quantum statistical approach to give the basic relations, for an extended 
description see Ref. [15]. We consider only strong interaction and neglect /3 equilibrium due 
to weak interaction. In addition to temperature, the total numbers of protons and neutrons as 
conserved quantities are fixed. In the grand canonical ensemble, the corresponding chemical 
potentials p„,p. p are introduced. The relation to the nucleon numbers is given by the equation of 
state 



where Q. denotes the system volume. The single nucleon state \p, cr, t) contains the momentum 
p and the isospin t — n,p. The summation over spin cr gives the factor 2. The Fermi function 
reads f T (u>) = [exp((w - p T )/T) + and A T (p, u>) is the spectral function of the nucleons. 
The spectral function A T (p, u>) is related to the self-energy "E T (p, cS) according to 



where E T (p) = p 2 /2m T (non-relativistic case). For small ImS r (//,w - ie), a peak arises at the 
quasiparticle energy Ef{p) - E T (p) + ReX T (p, Ej P (p)). The single-nucleon quasiparticle energy 
that depends on the properties of the nuclear medium is usually taken in the effective mass ap- 
proximation Ej V (p; T, n, Y p ) = AE^ E (T, n, Y p ) + p 2 /2m*(T, n, Y p ) where AEf E denotes the (rigid) 
self-energy shift and m* the effective mass of the nucleons. 

In warm nuclear matter, significant contributions to the spectral function are due to the imagi- 
nary part of the self-energy. Within a cluster decomposition, see lfl5ll . the A-nucleon propagators 
arise. The corresponding A-nucleon spectral function describes also the formation of bound 
states, that appear as peaks at Ea, v (P) = £U,v + P 2 /(2Am) in the low-density limit. The index v 
denotes the internal quantum state of the A-nucleon cluster, P its c.o.m. momentum. In general, 
the spin and internal excitation states have to be considered to specify the internal quantum num- 
bers. This way, the correct low-density limit of warm nuclear matter as a mixture of different 
components, the nuclei, is obtained. 

The A-nucleon spectral function is modified with increasing density. Similar to the quasi- 
particle peak in the single-nucleon spectral function, the peaks in the A-nucleon spectral func- 
tion are shifted to quasiparticle energies Ef (P) = Ea, v {P) + Re T,a, v (P, ET v (P))- It is clear in 
other many-particle systems, such as plasmas or condensed matter, that bound states behave like 
new components and can be treated in the same way as elementary particles (chemical picture). 
Notice that the introduction of the quasiparticle concept for the light clusters is related to the 
behavior of the corresponding spectral functions, irrespective of any particular approximation, 
and can be related to measurable properties. 

To evaluate the self-energy, perturbation theory and diagram techniques can be used. The 
cluster expansion of the self-energy [15] yields the EOS 




(1) 



A T (p, lj) = 



21m "L T (p, cl> - ie) 



(2) 



[cj - E T (p) - Re Z T (p, to)] 2 + [Im T, T (p, u> - ie)] 2 ' 




A,v,P 
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n„(T,p p ,fi„) = i^(A-Z)/ A , z [£^(P)], (3) 



where 

/ A ,zM = [exp(£( W - Zn p - (A - Z)p n )) - (- 1)V (4) 

is the Fermi or Bose distribution function that depends on the inverse temperature ft = l/iksT) 
and the chemical potentials p p ,p„ (instead of the isospin quantum number we use the charge 
number Z). The internal quantum number v denotes the excited states of the cluster A, Z, includ- 
ing the continuum states. In addition to the free nucleons A = 1, where Z = for n and Z = 1 for 
p, all higher clusters are included. 

We arrive at the NSE in the low-density limit where the quasiparticle energies E*f (P) can be 

replaced by the energies Ea, v {P) of free nucleons and clusters. The account of the contribution 
of the continuum (scattering states) in the sum over v leads for A = 1 , 2 to the Beth-Uhlenbeck 
formula for the second virial coefficient ]2[]. 

Quasiparticle shifts have to be considered in warm nuclear matter for densities above 1CT 4 
frrT 3 . Single-nucleon quasiparticle energy shifts Ef{p) have been extensively investigated. Var- 
ious approaches are used such as Skyrme, RMF, or DBHF approximations. These are summa- 
rized in Sec. [T]and will not be discussed here further. We elaborate expressions for the cluster 
quasiparticle energies E*f (P), A > 1. Starting from a Bethe-Salpeter equation for the A-nucleon 
propagator, the treatment of light clusters in warm nuclear matter is traced back to the solution 
of the in-medium Schrodinger equation 

2 // matter (l . . .A; 1' . . . A')^ v p(l' ■■■A') - Ef v (P)^ AvP {l ...A), (5) 

V.. A' 

where // matter ( 1 ... A; 1 ' ... A') is the instantaneous part of the in-medium Hamiltonian. In cluster 



mean-field approximation, explicit expressions are given in Refs. ll5U17tll8ll . 



3. Light nuclei in matter: The in-medium effective Schrodinger equation 

The few-body problem describing A < 4 nucleons in hot and dense matter can be related to 
an in-medium wave equation (Bethe-Salpeter equation) that is derived from many-particle ap- 
proaches, see Ref. 11511 . We consider only bound states and include A in the quantum state v, 
dropping the spin quantum number. For the light elements, v = d,t,h,a denotes the deuteron 
( 2 H), the triton ( 3 H), the helion ( 3 He), and the a particle ( 4 He). Considering uncorrected nu- 
cleons in the medium, the few-nucleon wave function and the corresponding eigenvalues follow 
from solving the in-medium Schrodinger equation 

[£f(l) + ... + ^ p (A)]<AvHl...A) 

+ Z Z [1 ~ - fiOWW, i'f) fl Wvpd' ■■■A') 

\'.. A' i<j kj=i,j 

= £f(F)<Ml...A). (6) 

For brevity, the single-nucleon quasiparticle energy Ej^(pi) is denoted as £'j P (l). The nucleon- 
nucleon interaction V(ij, i'f) becomes medium dependent due to the Pauli blocking prefactor 
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[1 - fi(i) - fi(f)]. The phase space occupation is described by a Fermi distribution function 
normalized to the total density of nucleons, 



1 n T 1 2nh 2 

= expt^D/r-^/n + i * T \kt 

in the low-density, non-degenerate limit (p T < 0). The chemical potential /i T is determined by the 
normalization condition Yipfiip) - n r, where t denotes isospin (neutron or proton), and 
has to be expressed in terms of these densities and the temperature. We restrict our calculations 
to the low-density region where the non-degenerate limit can be used, and replace the effective 
nucleon masses m* by the averaged free nucleon mass m so that h 2 /m = 41.46 MeV fm 2 . 

The in-medium Schrodinger equation © contains the effects of the medium in the single 
nucleon quasiparticle shift as well as in the Pauli blocking terms. Obviously, the bound state wave 
functions and energy eigenvalues as well as the scattering phase shifts depend on temperature and 
density. In particular, we obtain the cluster quasiparticle shifts 

Ef{P; T, n, Y p ) - E V (P) = AEf(P; T, n, Y p ) + AE^iP; T, n, Y p ). (8) 

Of special interest are the binding energies E* mi (P) = Ey P (P) - AE qp (P/A) that indicate the en- 
ergy difference between the bound state and the continuum of free (scattering) states at the same 
total momentum P. This binding energy determines the yield of the different nuclei according 
to Eq. OJ, where the summation over P is restricted to that region where bound states exist, i.e. 
Ef nd (P) < 0. 

The contribution of the single nucleon energy shift to the cluster self-energy shift AE^ E is 
easily calculated in the effective mass approximation, where the single-nucleon quasiparticle 
energy shift A£^ E (1) can be represented by the energy shift AE^ E and the effective mass [m*] -1 = 
[m T ] _1 + d 2 AE^ p (p)/dp 2 \ p= o. In the rigid shift approximation where m* = m, the energy shift 
A£f E cancels in the binding energy and can be absorbed in the chemical potential of the EOS 
Q. The effective mass correction is also easily calculated for given wave functions, see 1113111511 . 
The influence of the self-energy shifts on the binding energy is small and will be considered 
below in Sec. 14.31 

We consider here the Pauli blocking shift of the binding energies AEy* uh (P; T, n, Y p ) = 
Ef(P)-E v (P)-(A-Z)AEf i -ZAEf i . To evaluate it, we need the interaction potential V(ij, i'f). 
We choose a simple separable interaction potential with Gaussian form-factors 

(Pi-Pi? > 2 

V(p i ,pj;p' i ,p' j ) = Ay8p. +P]iP ' i+P '.e e *?* . (9) 

The two parameters A Y and y v can be fitted to reproduce the binding energy and the root mean 
square (rms) radii of the free nuclei as empirical input Ei. 

The solution of the few nucleon problem in the low-density limit can be found from varia- 
tional, Faddeev i2lll22ll . Green's-function Monte Carlo ll20ll . etc., approaches. We use a varia- 
tional approach with the Jastrow ansatz [ 15] for >J/ V p(l . . . A) that reproduces the exact solution 
for A = 2, 

1 1 — r P 

The prefactor Af„ is determined by normalization. 
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Table 1 : Light cluster wave function parameter values at zero density related to the Jastrow approach Eq. )10t 



V 


A v 




fly 


b v 




plan 


rmSy point 




[MeV fm 3 ] 


[fm- 1 ] 


[fm- 1 ] 


[fm- 1 ] 


[MeV] 


[MeV] 


[fm] 


d( z H) 


-1287.4 


1.474 


1.474 


0.2317 


-2.225 


10.338 


1.96 


f( 3 H) 


-1467.0 


1.153 


1.595 


0.567 


-8.482 


23.735 


1.68 


h ( 3 He) 


-1431.9 


1.153 


1.602 


0.5514 


-7.718 


23.021 


1.71 


a ( 4 He) 


-1272.9 


1.231 


2.151 


0.912 


-28.30 


51.575 


1.45 



Details are given in Ref. lfl5ll where Jacobian coordinates are introduced. The fit of the 
potential and the corresponding wave functions are given in Tab. Q] Due to the tensor force, the 
nucleon-nucleon interaction in the spin-triplet state is stronger than in the spin-singlet state. As 
detailed in Appendix A different orbitals for protons and neutrons will occur for A — 3, and 



different rms radii for protons and neutrons are considered there. The present analysis is based 
on averaged rms radii for t, h given in Tab.Q] 



4. Binding energies at finite densities 

4.1. Pauli shift in the low-density limit 

We focus on the Pauli blocking shift AEy ^(P; T, n, Y p ) that is responsible for the disap- 
pearance of bound states. It is calculated from Eq. (|6]l when the influence of the medium is 
accounted for only in the Pauli blocking terms. 

The solution is found in the low-density limit by perturbation theory, at arbitrary density 
by variational calculations. We use the Jastrow ansatz ( TTOb and determine optimal values for 
fly and b v . The evaluation of normalization, kinetic and potential energy is performed using 
Jacobi coordinates lfl5ll . Note that the in-medium Hamiltonian (O is not hermitian but can be 



transformed to a hermitian Hamiltonian, see Appendix C 

The evaluation of the Pauli blocking shift A£ , y ' aull (f > ; T, n, Y p ) is rather involved and time 
consuming. In the variational approach, we have to calculate multiple integrals and to search 
for a minimum. Similar to the single nucleon quasiparticle shift, where instead of the more 
fundamental DBHF calculations a simple fit within the RMF approach is convenient for further 
calculations, we will give simple expressions for the Pauli blocking shift that reproduce the 
numerical evaluations. 

In the low-density limit, the perturbation in Eq. (O is proportional to the neutron density 
(1 - Y p )n and the proton density Y p n. We discuss first the dependence on the asymmetry Y p . 
For the deuteron, Pauli blocking is determined by the sum of the neutron and proton distribution 
functions f„(pi) + f P (pi), see Eq. ([6j. It depends only on the total baryon density n in the non- 
degenerate limit. Therefore, we neglect the dependence on Y p for v - d. The same applies for 
4 He because neutron and proton orbitals are equally occupied for v - a. In the clusters with 
A — 3, however, neutrons and protons contribute differently to the internal structure, so that the 
shifts of 3 H and 3 He are sensitive to the asymmetry of nuclear matter. We find 

AE* m \P; n,T,Y p ) = f v (P;T,0) y v (Y p ) n + 0(n 2 ) , (11) 
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Table 2: Parameter values for the Pauli blocking shift A£^ allll (P; T, n, Y p ), Eq. 1121 , in the low-density limit 



V 


/v,0 


/v,l 


hi 


/v,3 


fv,4 




[MeV fm 5/2 ] 


[MeV fm 3 ] 


[MeV] 




[fm- 1 ] 


d( 2 H) 


388338 


6792.6 


22.52 


0.2223 


0.2317 


f( 3 H) 


159080 


20103.4 


11.987 


0.85465 


0.9772 


h ( 3 He) 


153051 


19505.9 


11.748 


0.84473 


0.9566 


a ( 4 He) 


352965 


36146.7 


17.074 


0.9865 


1.9021 



wi\hy d (Y p ) = y a {Yp) = 1, fortritony,^) = (| - jY p ), and for helion y h (Y p ) = (§ + §7 p ). For 
example, in comparison with helions ( 3 He), the tritons ( 3 H) show a stronger shift in neutron-rich 
matter because the neutrons in the cluster are stronger blocked than the protons. 

The functions f v {P; T, 0) can be calculated in first order perturbation theory using the unper- 
turbed wave functions of the free nuclei. Motivated by the exact solution for A = 2, we use the 
following fit for arbitrary v: 



f v (P; T,n) = f v ,i exp 



P 2 /h 2 



xlm < exp 



4(/ v y/ v 3)(l + 77/v,2) + «v» 



p/h 



1 2/ v ,4 

t 1 ' 2 P/h 



xerfc 



/v,3(l+/v,2/T) 



1/2 



1 -i- 



2/ Vj4 (i + r// V;2 ) 
p/h 

"2/ v ,4(i + r// y , 2 ) 



(12) 



where the term m v « can be neglected in the low-density limit considered in this subsection. 
In particular, at zero c.o.m. momentum P — we have 



/v(0; T,0) = 



(fy,2 + TW 



F( Xv ), x v = / y , 3 (i+/ v , 2 /r) 



1/2 



(13) 



where F(x) = 2x 2 (l — 7r 1 ^ 2 ;ce Jc2 erfc(x)). The complementary error function is defined as erfc(;t:) = 
1 - erf(jc) = 1 - 2tT 1/2 £ exp(-f 2 )df so that lim. MOO F(x) = 1 - 3/(2x 2 ) + . . ., lim A ^ F(x) = 
2x 2 (l - n l l 2 x + 2x 2 - n x ^ 2 x^ + ...). Expression ( [TBI is exact for the two-nucleon case, where (in 
units of MeV, fm) a d = y d = IMA, b d = f dA = (-E b d ind m/h 2 ) 1/2 = 0.2317, f da = h 2 a 2 d /(4m) = 
22.52, f di3 = 2 l/2 b d /a d = 0.2223, 



fd.O 



2\ 5 /2 



2 1 / 2 7T 3 / 2 fl 2 



and 



fd,i= - 



2\3/2 



2-(2 + l// 2 3 ) J F(^, 3 ) 



^ 2 2 3/2 / d ,3 



388338, 



= 6792.6 . 



1 -(1 + 1/(2/^) W43) 

The corresponding parameter values for the other light nuclei are given in Tab. [2] The intended 
accuracy to reproduce the solution of Eq. <JSJ is better than 1 % in the parameter region considered 
here. 
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Note that the dependence of the Pauli blocking shifts on the asymmetry parameter Y p be- 
comes more involved if we consider different orbitals for the neutrons and the protons in the 
clusters. In particular, this applies for triton ( 3 H) and helion ( 3 He) where the neutron and proton 



orbitals are different, see Appendix A 



4.2. Pauli shift at arbitrary subsaturation density 

The results for the Pauli shift obtained in the low density limit using perturbation theory can 
be extended to higher densities using a variational approach to solve the in-medium Schrodinger 
equation Eq. ©. We used the Jastrow ansatz ( fTOb for the class of wave functions and optimized 
the parameter values for a v , b Y . 

In the low-density limit, a linear dependence of the energy shifts on the nucleon density n 
follows from perturbation theory. For a more general dependence on the total nucleon density 
we consider the expression 



A£^ auli (P;n, T,Y P ) 



c y (P; T)il- exp 



MP; T,n) 

' c v (P;T) 



y v (Y p )n - d v (P; T, n)n" 



(14) 



This ansatz has been taken to ensure that the shift will not decrease with increasing density 
and to avoid spurious reappearance of bound states at high densities, after they are blocked out. 
The linear term f v (P; T, 0) is given by first order perturbation theory with respect to the density. 
The asymptotic term c v (P; T) should be larger than the cluster binding energy so that the bound 
state merges with the continuum of scattering states at a certain density and disappears. We 
restrict ourselves to the quadratic term with respect to the density dependence, a further density 
dependence in the prefactor d v (P; T, n) disappears at P — as discussed below. 

We have reduced the parameter dependences as much as possible to get the intended accuracy 
of few percent to reproduce the Pauli shifts also at higher densities. For zero momenta, P — 0, 
the temperature dependence of c v (0; T) and d v (Q; T, ri) is expressed as 



c v (0; T) = c vfi + 



Cy.\ 



d v (Q; T, n) = 



dy\ 



(T - dy^) 2 + c/ v ,3 ' 



The corresponding parameters are given in Tab. [3] 

The Pauli shift at finite momenta is fitted with c v (P; T) = c v (0; T) not depending on P, but 



(15) 



d v (P; T, ri) = d v (0; T, ri)e~ 



(16) 



An additional dependence on n, T is considered at finite values of P. Another additional de- 
pendence on n, T for finite momenta is introduced in f v (P; T, ri) where the dispersion relation 
becomes density dependent due to the parameter u v . Parameter values are given in Tab. [3] Start- 
ing from the fit to deuterons, the parameter values u v describing the dispersion were scaled to A 2 , 
see also |[l2ll where similar values have been considered. The optimalisation of the parameter 
values v v was performed to reproduce the calculated values for the Pauli blocking shift within 
few percent. 

4.3. Self-energy shift 

The quasiparticle shifts, Eq. ([§), contain besides the Pauli blocking also self-energy terms. 
The self-energy shifts can be easily evaluated if the effective mass approximation is taken for the 
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Table 3: Parameter values for the Pauli blocking shift A£y ™(0; T, n, Y p ), Eq. )14t , in units of MeV and fm (c Vj o, Cy,2>^v,2 
- [MeV]; c Vj] - [MeV 3 ]; c v , 3 ,d y , 3 - [MeV 2 ]; <f v ,i - [MeV 2 fm 6 ]; u v - fm; v„ - [MeV" 1 fm]) 



V 


d( z H) 


ten) 


h ( j He) 


a ( 4 He) 


Cy,0 


2.752 


11.556 


10.435 


150.71 


Cy,l 


32.032 


117.24 


176.78 


9772 


Cy,2 





3.7362 


3.5926 


2.0495 


Cy,3 


9.733 


4.8426 


5.8137 


2.1624 


d v ,i 


523757 


108762 


90996 


5391.2 


d v ,2 





9.3312 


10.72 


3.5099 


d v ,3 


15.273 


49.678 


47.919 


44.126 


u v 


11.23 


25.27 


25.27 


44.92 


Vy 


0.145 


0.284 


0.27 


0.433 



single nucleoli quasiparticle dispersion relation. Taking the unperturbed Jastrow wave function 
( [Tol l, the perturbative treatment gives the result 

(17, 

where the kinetic energy E^ m of the internal motion of the cluster v is given in Tab. [JJ It results 
as the averages of h 2 /mq 2 for A = 2, H 2 /m(q 2 + 3/4q^) for A = 3, and H 2 /m(q 2 + 3/Aq\ + 2/3qf) 
for A = 4, where g, denote the respective Jacobian momenta. To estimate the effect, an empirical 
expression m* jm — 1 - 0.17«/« sat can be taken. 

We denote the momentum P^° a (T, n, Y p ), where the bound state disappears, as Mott mo- 
mentum, see | Appendix B| At n > n^ ott (T, Y p ), the summation over the momentum to calculate 
the bound state contribution to the composition is restricted to the region P > P^ ott (T, n, Y p ). 
Below the Mott density, the influence of the self-energy shifts is small. If the nucleon density 
approaches the saturation density, the suppression of clusters is enhanced due to the self-energy 
shifts. 



5. Discussion and outlook 

The context of this work is the calculation of the composition of warm nuclear matter cover- 
ing a wide range of parameter values T, n, and Y p . Within a quantum statistical approach, exact 
results in the low-density limit as well as quasiparticle approaches near the saturation density are 
reproduced. A generalized Beth-Uhlenbeck approach [7] can be used that leads to NSE and the 
virial EOS at low densities, but includes medium effects that become relevant for baryon number 
densities above ~ 10~ 4 frrT 3 . 

In the present work, we focus on the light elements d, t, h, a immersed in warm and dense 
matter. The quasiparticle description is an important prerequisite for their proper treatment in 
connection with thermodynamic or transport properties. The parametrization given here can 
be improved in different ways. The solution of the few-body problem accounting for medium 
effects, Eq. ©, can be performed using more sophisticated methods instead of the variational 
approach. The class of functions that represent the quasiparticle shifts as function of temperature, 
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density, and asymmetry parameter, can be extended and better adjusted to the solution of the few- 
body problem for the light elements. 

We restrict us to the light elements A < 4. An extension to heavier elements is possible but 
needs further considerations, see Ref. [8]. Since the following elements up to carbon are weakly 
bound, only small yields are expected. In particular, at low temperatures the formation of heavier 
clusters is of interest. Then, the proper treatment of Coulomb interactions and inhomogeneous 
solutions is inevitable As well known, in the region of thermodynamical instability 

the Coulomb interaction is responsible for structure formation representing both phases. 

Another issue is the treatment of excited states, in particular the continuum of scattering 
states, that are also included in the EOS (0. The contribution of scattering states has to be taken 
into account to reproduce the correct low-density limit of the second virial coefficient. For the 
in-medium treatment of scattering states in the two-nucleon case, as well as the evaluation of the 
second virial coefficient, see Qa LZl l2D • In particular, the deuteron fraction is reduced if scattering 
states are considered, because the binding energy is comparable to the temperature. Sharply 
peaked structures related to the Mott effect are washed out because of the Levinson theorem, see 
I2D. For nonequilibrium processes like HIC, the pole structure of the cluster spectral function 
may be relevant so that the cluster yields are not determined by the contribution of scattering 
states to the second virial coefficient. 

The extension of the low-density results for the quasiparticle shifts to the region of the satura- 
tion density needs further discussion. The expressions for the self-energy and the Pauli blocking 
are obtained in first order with respect to the distribution function /i (1) in the medium part of the 
Hamiltonian. Higher orders will arise, e.g., from the free A-particle propagator. Also, the peaks 
in the spectral function that characterize the quasiparticle excitations are broadened at increasing 
density what reflects the damping of the quasiparticles due to collisions. 

Moreover, the approximation of the uncorrected medium used in the in-medium wave equa- 
tion (|6]l can be improved considering the cluster mean-field approximation Ia, ll5lll7l1 . In partic- 
ular, at low temperatures and low densities a cluster become dominant, and the medium effects 
of a matter are produced by the surrounding a clusters. Then, the a quasiparticle energy shift is 
determined by the effective a - a interaction. Finally, degeneracy and the formation of quantum 
condensates (pairing, quartetting) has to be considered in the low temperature region, when the 
single-nucleon chemical potentials approach the lowest bosonic (A = 2,4) binding energy. 

In conclusion, the parametrization of the quasiparticle cluster energies provides us with a 
tool to extend low-density approaches such as NSE to higher densities where mean-field effects 
become relevant. This is an indispensable ingredient to calculate the composition of nuclear 
matter in astrophysical applications as well as in HIC. The contributions to the quasiparticle 
energies that are linear in the density, expressed by the functions f v (P; T), Eq. ( TBI , are well 
described already in perturbation theory. At higher densities, exploratory results have been given, 
in particular with respect to the disappearance of bound states with given c.o.m. momentum. The 
detailed treatment of correlations formed in the nucleonic many-body system at higher densities, 
in particular the evaluation of the spectral function, remains a challenging problem. 

Appendix A. Binding energies and rms radii of light nuclei 

We start out from the empirical values for the binding energies and charge rms radii for 
the light elements see Table IA.4I The experimental charge radii have been converted 

to point proton rms radii by removing the proton and neutron (r 2 ) of 0.743 and -0.116 fm 2 , 
respectively. 
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Table A.4: Empirical data of light nuclei GjJUi 





binding 


mass 


spin 


rms-radius 


rms-radius 




energy 






(charge) 


(point, proton) 




[MeV] 


[MeV/c 2 ] 




[fm] 


[fm] 


n 





939.565 


1/2 


0.34 





pCu) 





938.783 


1/2 


0.87 





d( 2 H) 


-2.225 


1876.12 


1 


2.12 


1.96 


f( 3 H) 


-8.482 


2809.43 


1/2 


1.76 


1.58 


h ( 3 He) 


-7.718 


2809.41 


1/2 


1.93 


1.76 


a ( 4 He) 


-28.30 


3728.40 





1.68 


1.48 



Because of equal numbers, the neutron and proton orbitals are approximately identical in d 
and a. The situation is more complex for t and h where the numbers of protons and neutrons 
are different. We give some exploratory calculations to estimate the different orbitals of neutrons 
and protons in t and h, in particular the different rms radii. 

We deduce characteristics of the nucleon wave functions from experimental data that are 
available for the binding energies E,,Ei,, as well as the proton rms radii rms,j, and rmsft iP . 
Coulomb interactions and tensor forces are responsible for the differences between t and h. As 
mirror nuclei, the difference in the binding energy is caused by the Coulomb interaction that leads 
to a repulsion between the protons in 3 He so that the binding energy is reduced if comparing with 
3 H. We also expect that the rms values are slightly larger in h because it is not so strongly bound 
as t. However, the observed large difference in the rms radii is due to the tensor force of the 
nucleon-nucleon interaction that is related to spin orientation. 

The occurrence of the deuteron as a bound state in the p—n spin triplet channel, in contrast to 
p — porn- n where a bound state is absent, is due to the tensor force that makes the interaction 
between nucleons in the spin-triplet state stronger than between nucleons in the spin-singlet state. 
In the ground state of 3 H, there are two neutrons with opposite spin orientation, whereas a proton 
is found, say, with spin up. The interaction of this proton with the neutrons contains contributions 
from the tensor force (with the spin-up neutron). The spin-up neutron feels also the tensor force 
(with the proton), whereas the spin-down neutron sees only opposite spins so that no tensor force 
acts. On average, the neutron is weaker bound in 3 H than the proton so that its rms radius is 
larger, rms, „ > rms, tP . Vice versa, in 3 He the proton rms radius is expected to be larger than the 
neutron one, rmsjj, > rms^. 

To separate both, the tensor and Coulomb effects, we assume first that Coulomb effects can 
be discarded. Then we have rms',„ = rms/ ljP and rms'/, „ = rms fjP . In the average, the point 
nucleon rms radius of the three-nucleon bound state is rms'3 aver = (rms fiP + 2 rms/, p )/3 = 1.70 
fm. Correspondingly, we consider an averaged binding energy E^ wm = {E, + 2£ , /,)/3 = 7.973 
MeV. This averaged three-nucleon bound state serves only for guidance to estimate the effects 
and does not provide high precision, because it is contaminated by Coulomb effects. For this 
average nucleus, we obtain the parametrization A3 - 1443.66,73 = 1.153, 03 = 1.5998, and 
b 3 = 0.5567. 

For the further analysis, we assume that the range of the interaction y is nearly constant. To 
study the influence of the Coulomb interactions we change the strength A of the interaction and 



11 



Table A. 5: Light cluster wave function parameter at zero density from the Jastrow approach, A = 3 



V 


Ay 


Jv 


fly 


b v 


Ey 


rms v point 




[MeV fm 3 ] 


[fm- 1 ] 


[fm- 1 ] 


[fm- 1 ] 


[MeV] 


[fm] 




-1600.5 


1.153 


1.571 


0.6237 


-11.525 


1.59 


t ( 3 H)n 


-1413.2 


1.153 


1.607 


0.5429 


-7.320 


1.73 


h ( 3 He)p 


-1381.5 


1.153 


1.615 


0.5281 


-6.645 


1.76 


h ( 3 He)n 


-1553.6 


1.153 


1.578 


0.6041 


-10.441 


1.62 



evaluate the change of binding energy and of the rms radius. This gives a relation between the 
binding energy and the rms radius. In particular, for £3 = E, we find rms, = 1.68 fm, and for 
Ej = E/, we find rms/, =1.71 fm. From this, we conclude that the weaker bound nucleus (h) is 
blown up by rms/, - rms, = 0.03 fm. We determine the potential parameters A,y that reproduce 
the values £,,rms,, and Zs/^rms/,, see TableQ] 

Note that we can specify the neutron point rms radius, assuming that the orbits in h are ex- 
panded by 0.03 fm compared to t. Using the measured rms radii for the protons, we estimate 
values for the neutron point rms radius rms^^rms/,,,, shown in the Table |A31 Finally, we de- 
termined the proton and neutron wave functions assuming that the range of interaction remains 
fixed and the strength of interaction is fitted to the respective rms radii. Results are given in 
Table I A. 51 Using the corresponding parameter values, we can calculate the contributions of the 
different orbitals to the Pauli blocking shift of the triton t and helion h. We will give the results 
in a forthcoming work and restrict us here to the shifts derived from the averaged values given in 
Tabled] 



Appendix B. Pauli blocking shift of the deuteron 

In the case A — 2, we introduce Jacobian wave numbers q\ — (p2 - pi)/2h,hq2 — P — 
p\+P2- The exact solution of Eq. (0 gives the momentum dependent in-medium binding energy 
E* md (P) = Ef{P) - P 2 /(4m) = E d + AE p / ali (P) according to (K = P/K) 

— r dq\ az\—— — — e y <i 

4tt 2 J,, * J_i E^ A (P)-W-q\jm 



(B.l) 



with effective chemical potentials (containing the self-energy shifts) determined by 

1 r°° 1 

dkk 2 — =n T . (B.2) 

Numerical solution gives the binding energy E^ lnd (P; T, n, Y p ). 

In the non-degenerate case, i.e. at low densities and not too low temperatures, where the 
effective chemical potentials are negative, integrations can be performed analytically. Then, it is 
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0.005 0.01 0.015 0.02 

baryon density [fm ] 



Figure B.l: Shift of the deuteron binding energy £* M (P = 0) = E d + AE p / uh (Q), Eq. {8), as function of the baryon 
density n at different temperatures T, Y p = 0.5. The numerical solution of Eq. )B.U is compared with the fit formula 
(thin lines), Eq. 1141 . and the non-degenerate approximation (dotted lines). 



easily shown that only the total baryon density n — n„ + n p enters so that the solution does not 
depend on the asymmetry Y p . The typical behavior of the binding energy is shown for symmetric 
matter (Y p = 0.5) in Figs. IB.ll and lB.2l For comparison, the solution replacing the Fermi function 
by the Boltzmann function is also shown. At the temperatures considered, the differences are 
small. Only at very low temperatures where Bose-Einstein condensation is possible, significant 
differences are expected. 

Depending on temperature, the absolute value of the binding energy at zero momentum 
£bmd(Q. n ^ j s decreasing with increasing density and merges with the continuum of scat- 
tering states. The analytical fit reproduces this behavior with exception of the region where the 
bound state merges with the continuum. We denote this density as Mott density n¥ ott (7', Y p ). 
Note that this disappearance of a bound state does not introduce discontinuities in the EOS, 
provided that the contribution of scattering phase shifts is also taken into account |7|]. 

To give a simple approximation for the Mott density we perform a fit (in units MeV, fm~ 3 ), 

nf oa (T, Y p ) « 0.00057646 + 0.00065443 T + 0.000012491 T 2 . (B.3) 

For baryon densities higher than the Mott density, bound states can be formed only for c.o.m. 
momenta P larger than the Mott momentum hK^ lott (T, n, Y p ). Using the solution of Eq. ( IB .It at 
arbitrary degeneracy, a fit is given as 

\K Mon (T nY)f~ 4 - 5185 - Q - 16164r + ™ 2r2 
L d 1 1 ~ 2(1.32- 0.02782 T) 

( ( 4.5185 - 0.16164 T + 0.0056582 T 2 ) 2 1000(» - nf*<(T, Y p ) )) l/2 

4(1.32 - 0.02782 T) 2 + 1.32- 0.02782 T 
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deuteron c.o.m. wave number [fm ] 



Figure B.2: Shift of the deuteron binding energy E^ ni (hK) = Ej + AEj m] \hK), Eq. {8j, as function of the deuteron 
c.o.m. momentum P = hK at T = 10 MeV for different baryon densities n, Y p = 0.5. The numerical solution of Eq. IB.lt 
is compared with the fit formula (thin lines), Eq. 1141 . The dotted line is the fit according Eq. jB.6t . 
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(B.4) 

The behavior of E^ md (P; T, n, Y p ) as function of the c.o.m. momentum P is shown in Fig. IB. 21 
for different densities at fixed temperature T = 10 MeV and Y p = 0.5. With increasing P, the 
overlap of the wave function with the occupied phase space becomes smaller so that the Pauli 
blocking mechanism becomes less efficient, and the shift is decreasing. A Gaussian fit as used in 
Refs. lfl2l[l5ll is reasonable at small densities. The merger with the continuum leads to deviations 
so that, in principle, higher orders of P 2 can be considered in the dispersion relation. 

To give a more general fit that reproduces the exact solution within a few percent within the 
entire P space, we take 



AE p / u,i (P; T, n, Y p ) = A£™"(0; T, n, Y p ) 



Pauli/ 



e -gaAT,n)P 2 /tl 2 

l+g b , d (T,n)P 2 /h 2, 



for n < n™ ott (T, Y p \ and for n > n™ ott (T, Y p ): 



AE p / ah (P;T,n,Y p ) = E d 



l+gcAT,n)(K-Kf° a (T,n,Y p )) 



with the quantities 

gaA( T < ") 

gb,d(T> ") = 



(0.279145 - 0.0099364 T + 0.00017819 T 2 ) 
+(0.41393 + 0.0069374 T)n- 1.65 n 2 , 



8mT 



1 - 



2 fdlfd,3 

3T 



1 



F(f d ,3(i+fdi/Ty 2 ) 

: [l + (25.797 - 10.456 T + 0.29668 T 2 ) n] . 



- 1 



gc,d( T > ») = 



1 - 



1 



- 1 



- ga,d( T > «) 



- ga,d( T > «) 



8mT [' 3T [F{f d ,i{\ + f d , 2 ITfl 2 ) 
x [l - (41.703 - 2.6883 T + 0.18684 T 2 ) n] 
+(1.0613 - 0.083155 T + 0.0018922 T 2 )Kf° tt (T,n, Y p ) . 



(B.5) 



(B.6) 



(B.7) 



Large deviations of the fit and the exact solution arise for small c.o.m. momenta near the Mott 
density n™ ott (T, Y p ), see Fig. ~ 



Appendix C. Symmetrization of the in-medium Hamiltonian 

The in-medium Hamiltonian © is not hermitian, H mMer (l ...A,Y .. .A') + [H matter (l' . . . A', 1 . 
see also rf23ll . In the case A = 2 it can be transformed to a hermitian Hamiltonian, [I - f(l) — 
/(2)]w(l,2)w(l'2') [1 -/(l)-/(2)] 1/2 w(l,2)w(l'2')[l -/(l')-/(2')] 1/2 transforming also 
the wave functions \l/dp{\2) — > [1 - /(l) - /(2)] _1 ^ 2 i/rj/)(12). In the nondegenerate case where 
/(l) <s; 1, we can take [1 - /(l)/2 - /(2)/2]w(l, 2)w(l'2')[l - /(l')/2 - /(2')/2] for the sym- 
metrised Hamiltonian. In the degenerate case, new effects such as pairing will occur. In contrast 
to the Pauli blocking term [1 — /(1)][1 - /(2)] used in the Brueckner theory, where hole-hole 
contributions are neglected, the full Pauli blocking term [1 - /(l) - /(2)] leads to the Gor'kov 
equation in the BCS theory of superfluidity. 

15 



References 



[1] A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle systems (McGraw-Hill, New York, 1971). 
[2] J. M. Lattimer and F. D. Swesty, Nucl. Phys. A 535, 331 (1991). 

[3] H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi, Progr. Theor. Phys. 100, 1013 (1998); Nucl. Phys. A637 435 
(1998). 

[4] T. Klahn et al, Phys. Rev. C 74, 035802 (2006). 

[5] A.S. Botvina and I.N. Mishustin, ArXiv e-prints. larXiv:0811.2593"l (2008). 

[6] G. Ropke, L. Miinchow, and H. Schulz, Nucl. Phys. A 379, 536 (1982); Phys. Lett. B 110, 21 (1982); G. Ropke, 

M. Schmidt, L. Miinchow, and H. Schulz, Nucl. Phys. A 399, 587 (1983). 
[7] M. Schmidt et al, Ann. Phys. 202, 57 (1990). 

[8] G. Ropke, M. Schmidt, and H. Schulz, Nucl. Phys. A424, 594 (1984); G. Ropke, Clusters in Nuclear Matter, in: 
Condensed Matter Theories, vol. 16, Ed. S. Hernandez and J. W. Clark (Nova Sciences Publ., New York 2001); 
Phys. Lett. B 185, 281 (1987). 

[9] C. J. Horowitz and A. Schwenk, Nucl. Phys. A 776, 55 (2006). 

[10] G. Shen, C. J. Horowitz, and S. Teige, Phys. Rev. C 82, 045802 (2010); ArXiv e-prints. larXiv:1101.3715l (2011). 

[11] M. Hempel and J. Schaffner-Bielich, Nucl. Phys. A 837, 210 (2010). 

[12] S. Typel et al, Phys. Rev. C 81, 015803 (2010). 

[13] K. Sumiyoshi and G. Ropke, Phys. Rev. C 77, 055804 (2008). 

[14] J. Natowitz et al, Phys. Rev. Lett. 104, 202501 (2010). 

[15] G. Ropke, Phys. Rev. C 79, 014002 (2009). 

[16] S. S. Avancini, L. Brito, J. R. Marinelli, D. P. Menezes, M. M. W. de Moraes, C. Providencia, and A. M. Santos, 
Phys. Rev. C 79, 035804 (2009); C. J. Horowitz and D. K. Berry, Phys. Rev. C 78, 035806 (2008). 

[17] G. Ropke, T. Seifert, H. Stolz, andR. Zimmermann, Phys. Stat. Sol. (b) 100, 215 (1980); H. Stolz, R. Zimmermann, 
and G. Ropke, Phys. Stat. Sol. (b) 105, 585 (1981). 

[18] J. Dukelsky, G. Ropke, and P. Schuck, Nucl. Phys. A 628, 17 (1998). 

[19] H. de Vries, C. W. de Jager, and C. de Vries, Atomic Data and Nuclear Data Tables 36, 495 (1987). 

[20] S. C. Pieper, R. B. Wiringa, Ann. Rev. Nucl. Part. Sci. 51, 53 (2001). 

[21] M. Beyer, W. Schadow, C. Kuhrts, and G. Ropke, Phys. Rev. C 60, 034004 (1999). 

[22] A. Sedrakian and J. W. Clark, Phys. Rev. C 73, 035803 (2006); A. Sedrakian and G. Ropke, Ann. Phys. (NY) 266, 
524 (1998). 

[23] P. Danielewicz, Ann. Phys. 197, 154 (1990); D. Kremp, M. Schlanges, and W.-D. Kraeft, Quantum Statistics of 
Nonideal Plasmas, p. 195 (Springer, Berlin, 2005). 



16 



